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^vQ , We present a mean-field theoretical study on the effect of a single non-magnetic impurity in quasi- 

one dimensional unconventional density wave. The local scattering potential is treated within the 
self-consistent T-matrix approximation. The local density of states around the impurity shows the 

1 \ presence of resonant states in the vicinity of the Fermi level, much the same way as in d-density 

waves or unconventional superconductors. The assumption for different forward and backscattering, 
characteristic to quasi-one dimensional systems in general, leads to a resonance state that is double 
peaked in the pseudogap. The Friedel oscillations around the impurity are also explored in great 
detail, both within and beyond the density wave coherence length fo- Beyond £o we find power 
i— H ' law behavior as opposed to the exponential decay of conventional density wave. The entropy and 

specific heat contribution of the impurity are also calculated for arbitrary scattering strengths. 

H ■ 

^ . PACS numbers: 71.45.Lr, 75.30.Fv, 71.27. +a 
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The effects of a single magnetic and non-magnetic impurity on the properties of high temperature superconductors 
(HTSC) have attracted considerable attention in the past few years both from the theoretical?^' 3 -'^' 5 -' 6 - and more 
recently from the experimental side£' 7 -'£' 9 - This enthusiasm from the theory is mainly due to the fact, that electronic 
scattering from a localized impurity produces qualitatively different electronic structure around the impurity in the 
d-wave superconductor phased (dSC) and in the alternative models proposed for the pseudogap phase of high- 
T c cuprates. These models include for example the d-wave density wav o 1 ' 3 ' 10 ' 11 (dDW), and the superconducting 
. phase-fluctuation scenario^ 2 (PF). As the local density of states (LDOS) is proportional to the differential tunneling 
^C) ' conductance, scanning tunneling microscopy (STM) measurements are in principle able to distinguish between the 
f*"** , different available scenarios due to their spatial and energy resolution capability. Indeed, recent STM spectroscopy 
t-H ■ performed on I^S^CaCr^Os+a HTSC below the superconducting transition temperature T c confirmed the existence 
of the impurity induced quasiparticle subgap resonance close to the Fermi energy, anticipated in a dSCi&S 

The understanding of the microscopic origin of the pseudogap (PG) phase of HTSC represents a formidable challenge 
for both theory and experiment. Many experiments have provided evidence for unusual normal state behavior in 
underdoped cuprates. One of the most striking is definitely the observed gaplike feature of the quasiparticle density 
of states (DOS) above T c , with depleted electronic states around the Fermi energy^ It was shown lately, that 
independently of the origin of the PG state, such a DOS is alone sufficient to produce a resonant impurity state close 
to the Fermi level in the presence of a local non-magnetic impurity^ The development of this single quasiparticle 
£h ', peak in the pseudogap violates particle-hole symmetry. Similar conclusions were reached in the dDWi 3 . and PF 3 
O ■ models too, although in the latter case the subgap resonance was found to be double peaked. 

The study of electronic scattering from a localized perturbation however is not limited to HTSC physics. In fact, 
this issue was on agenda already in the late sixties, in the context of conventional BCS superconductivity 15 : 16 ' 17 It was 
found that in s-wave superconductor both magnetic^ 5 - and resonant non-magnetioi^ impurities produce true bound 
states inside the energy gap. In addition to that, conventional charge density waves (CDW) were also investigated 
from this aspect ) 18 ' 19 : 20 ' 21 and intragap bound states were obtained here as well. 

In this theoretical study, our aim is to explore the effects of a single non-magnetic impurity in quasi-one dimensional 
unconventional density wave (UDW). As the scattering is spin independent, our theory and conclusions are applicable 
to both unconventional charge density wave (UCDW) and unconventional spin density wave (USDW) . An UDW is a 
density wave with momentum dependent gap A(k) t 22 ' 23 the average value of which over the Fermi surface vanishes. 
Thus these systems lack spatial modulation of either charge or spin density that could be observed in principle by 
conventional means (x-ray or NMR). Though this sort of hidden order makes experimental detection more complicated, 
still a great deal of effort has been done for example in organic conductor a -(BEDT-TTF) 2 KHg(SCN) 4 , 2 ii 2 £' 2 £i 2 l' 2 £i2 9 . 
The low temperature phase of this salt has been explained rather well by a quasi-one dimensional UCDW, see Ref. [22j 
and references therein. Potential scattering from impurities with finite concentration has been investigated in UDW in 
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the weak scattering Born limitjsP. in the unitary limit^ and more recently in an extended scheme valid for arbitrary 
scattering amplitudes 32 . To the best of our knowledge however, the single impurity problem in UDW has not been 
addressed so far. Therefore we study this issue here. We explore the electronic structure and Friedel oscillations 
around the impurity. Besides, the specific heat contribution is calculated as well for arbitrary scattering strength. 

The paper is organized as follows: in Section 2 we present the formalism of the T-matrix approach. In Sections 
3-5 the formalism is applied to quasi-one dimensional unconventional density wave in order to obtain numerical and 
analytical results for the local density of states around the impurity, the induced Friedel oscillations, and for the 
thermodynamics, respectively. Finally, Section 6 is devoted to our summary and conclusion. In Appendix A further 
useful results are collected, those we make use of throughout the article. 

II. FORMALISM 

We begin with the mean-field Hamiltonian describing pure quasi-one dimensional density wave o 23 ' 33 

H o=^2 £,( k )( a l,a a K<? - 4-Q, CT a k-Q, CT ) + A CT (k)4 CT a k _ Qj(T + A;(k)aj,._ Q)0 .Ok )ff , (1) 

where a k and ak, CT are, respectively, the creation and annihilation operators for an electron in a single band with 
momentum k and spin a. The prime on the summation means that it is restricted to the reduced Brillouin zone as 
\k x — < k c , where k c is the momentum cutoff. The best nesting vector is given by Q = (2kF,ir/b,ir/c), with Uf 
being the Fermi wave number. Our system is based on an orthorhombic lattice with lattice constants a, b, c towards 
directions x, y, z. The system is highly anisotropic, the kinetic energy spectrum of the particles linearized around the 
Fermi surface reads as 



£(k) = vp(k x — hp) — 2tbCos(bk y ) — 2i c cos(cfc z ), (2) 

and the quasi-one dimensional direction is the x axis. For (U)CDW, the density wave order parameter A cr (k) is even, 
while in (U)SDW it is an odd function of spin index. In the case of conventional density waves the order parameter is 
constant on the Fermi surfaced As opposed to this, for unconventional condensates it depends on the perpendicular 
momentum and has different values at different points on the Fermi surface. The precise k-dependence is determined 
by the matrix element of the electron-electron interaction through the gap equation^ 

In the presence of a local impurity, the order parameter is no longer homogeneous in real space. A self-consistent 
calculation with the interaction included results in a spatially varying A(r)^ where the deviation from the bulk value, 
arising from pair-breaking effect, is confined to the vicinity of the impurity itself. The Zn substitution in cuprates 
is a good example of this. However, the phenomena discussed here are not expected to be altered significantly by 
introducing position dependent corrections^ Thus, we completely ignore this effect and stick with the homogeneous 
solution, the momentum dependence of which is chosen as A(k) = Ae l( ^ sin(6fc a ), 23 - where the phase <j> is unrestricted 
due to incommensurability. 

The interaction of electrons with the single non-magnetic impurity placed at the origin is described by the Hamil- 
tonian 

ffjr'/ ak +q ,. V (U(0) U(Q)\ ( a Ka \ 

Hl - V ^ Uk+q-q.J \U(Q) 17(0) ) ^k-q, J ■ {6) 

k,q,CT 

Here the summation over momentum transfer q is restricted to small values, V is the sample volume and we have 
neglected the small q-dependence of the matrix elements (i.e. the Fourier components of the impurity potential). 
This is because we are dealing with quasi-one dimensional density wave systems, whose Fermi surface consists of two 
almost parallel sheets at the points ±/cf- Consequently, two relevant scattering amplitudes can be distinguished: the 
forward U(0) and the backward scattering parameter U(Q), respectively. We believe that we can capture the essence 
of physics with this approximation. This is indeed the case, at least as long as only static quantities, single particle 
properties and thermodynamics is concerned. Investigating the effect of pinning, and dynamic properties such as the 
sliding density wave however, requires the inclusion of the small q contribution as wel L 34 ' 35 

It is worthwhile to emphasize at this point, that the allowance for different forward and backscattering in Eq. (J3]) 
constitutes a more realistic impurity physics than the usage of a somewhat artificial point-like scalar impurity (U (r) = 
US(r)), with U(Q) = U(Q) = U. Nevertheless, this restricted impurity model is widely used in literature concerning 
the single impurity problem in conventional charge density waves , 18 ' 19 unconventional superconductors with d-wave 
pairing symmetry) 4 ^ and more recently in the context of the pseudogap phase of HTSC ^ 2 ' 3 ' 14 We will see shortly 
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however, that the generalized approach with U(0) ^ U(Q) will reveal the true double peaked nature of the subgap 
resonance in UDW. These two virtually bound states with finite lifetime correspond to the infinitely sharp, and 
therefore well-defined bound states found in fully gapped CDW.— We suspect, that similarly to quasi-one dimensional 
density waves, analogous results for the counterpart resonance peak could be obtained in quasi-two dimensional d- 
density wave. 

The main part of the paper is devoted to calculating the effect of impurity on the spectral function A(r, uS) (in 
other words the local density of states), and the total electron density n(r). These quantities are directly related to 
the single particle Green's function. In view of the change in thermodynamic potential SQ, caused by the interaction, 
we will also determine later on the specific heat contribution of dilute impurities at arbitrary scattering strengths. As 
the interaction Hamiltonian in Eq. ([3]) is quadratic in fermion operators, SQ is thus related to the Green's function 
as well. According to all these, our first goal is to calculate the single particle propagator dressed by the scatterer to 
infinite order. In view of this, the rest can be obtained in a relatively straightforward manner. Hereafter we drop spin 
indices since they are irrelevant for our discussion on the effect of potential scattering, and for the sake of simplicity 
consider spinless fermions but our conclusions apply to both unconventional spin and charge density wave. 

In the usual way the electron field operator can be split into left (L) and right (R) moving parts as \&(r) = 
*i(r) + ty R (r)^2. where 

*«W = -7f E' elkrflk ' a = R (= L (= ( 4 ) 

v k 

and the prime on the summation means the constraint \k x — akp\ < k c in accordance with Eq. (TT]). Choosing a sharp 
cutoff is somewhat artificial and obvio usly model dependent. Nevertheless, we shall follow the treatment developed 
for conventional density wave in Ref. [18| . and the cutoff will be taken to infinity whenever it does not affect the 
essence of physics. In this respect, most quantities we evaluate in the paper are insensitive to the applied cutoff. 
The definition of Green's function including the effect of impurity is 

G a p{v,v'-r) = -<T T [* c (r,r)*J(r')])H, (5) 

where H refers to the total Hamiltonian Hq +H\. Applying standard equation-of- motion technique, Dyson's equation 
can be given in a matrix form for the Matsubara components as^ 

G Q/3 (r, r'; iw n ) = G° a0 (r, r'; iuj n ) + G° 7 (r, 0; iuj n )t lS Gsp{0, r'; iu n ), (6) 

where t 7 s — U(0)S-yg + £^(Q)&y,-<5> the zero superscript stands for the bare propagator, and summation is applied 
over indices occurring twice. Now the spectral function, the particle density and the change in the grand canonical 
potential are explicitly given as follows 

A(r, uj) = - -Im V" G a/3 (r , r ; iu n — >• w + iO) , (7) 

7T ' — ' 

a/3 

n(r) = ^G Q/3 (r,r;r = -0), (8) 

a/3 

and finally, SD, is calculated from the coupling constant integral as 

Sn= [ d\ lim {U{0)[G rr (0,0;t) x + G ll (0,0;t) x } + U{Q)[G rl (0,0-t) x + G lr {0,0-,t) x }}. (9) 
Jo T ^-° 

In the following sections we give detailed description of the quantities above in UDW, and will compare them to 
those of a fully gapped conventional DWi 18 ' 19 i 20 ' 21 Besides, emphasis will be put also on the fact, that the local density 
of states around the scalar impurity in UDW, in particular its intragap structure, greatly resembles that of the dDW 
phase obtained in Refs. [3I3I3]. 

III. LOCAL DENSITY OF STATES 

In this section we present the results for the spectral function making use of Eq. {?])• Dyson's equation for the 
Green's function can be readily solved as 

G Q/3 (r, r'; »w w ) = G^(r, r'; iu n ) + G° 7 (r, 0; iW „)T 75 (^„)G^(0, r'; iu n ), (10) 
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FIG. 1: The dimensionless local density of states A(r,u)/No is shown in UDW versus energy and position measured from the 
scatterer for m = (top left), m = ±1 (top right), m = ±2 (bottom left) and m — ±3 (bottom right). For the numerics 
NoU(0) = 0.2, NoU(Q) — 0.7, vfIzf/A — Virkp/tb — 10 and <f> — ir/3 were applied. The black lines on the top left figure were 
calculated from kpx = nnvFkF / + ^f^f) with n = 5, 6, 7 and 8 (from bottom to top). 

FIG. 2: Left panel: the local density of states in UDW is plotted on the m = chain, enlarged around the Fermi energy to 
better show the resonant states in the pseudogap. The other parameters are the same as on Fig. [Tj Right panel: the same 
quantity for fixed distances from the scatterer, kpx = 0.9 (dashed), 1.8 (dashed-dotted) and 4 (solid). 

where 

T l5 {iw n ) = {tll-G^O.O;^]- 1 }^ (11) 

is the T-matrix. Equation (fT0|) provides us the general dependence of G a p(r, r'; iui n ) on r and r', in the followings 
however, according to Eqs. ([30 we consider only diagonal components (in real space). Furthermore, in a quasi-one 
dimensional UDW besides the chain direction x, in which the model is continuous, perpendicular spatial dimensions 
are present as well offering the possibility for momentum dependent gap . 22 ' 23 This is to be contrasted with the strictly 
one-dimensional nature of normal CDW. 36 As the specific k-dependence of the density wave order parameter was 
chosen to be A(k) ~ sm(bk y ), the relevant perpendicular direction is y, and the model remains discrete in this 
variable. Consequently we take r = (x, mb, 0), where b is the corresponding lattice constant and m is an integer 
indexing parallel chains. 

It is clear from the structure of the T-matrix, that while analytically continuing to real frequencies as prescribed 
in Eq. I|7|). poles can emerge at certain well defined energies that describe bound states localized at the impurity site. 
This phenomenon is well understood in a normal metal for quite a long time^i There, the binding energy of these 
states are either above the upper or below the lower edge of the band. Bound states can also show up in conventional 
density wave, though the situation in this case is more complex because additional states can develop in the gap 
as wel h 18 i 19 ' 20 i 21 We will see in a little while that in unconventional density wave, though being gapless^ localized 
intragap states develop near the Fermi energy (see the region \u>\ < A in Figs. [I] and [5]). Albeit these states closely 
resemble to the corresponding results in fully gapped charge or spin density wave, their widths in energy are not 
infinitely sharp and therefore the identification as true bound states is not adequate. Rather we call them virtual 
or resonance states. In normal metal, this sort of accumulation of states in a broader energy range is referred to 
as a virtual state as wel\M- On the other hand, the notion of a resonance state is widely used in the terminology 
of unconventional (d-wave, etc.) superconductivity; 4 ^ and more recently in the context of the pseudogap phase of 
high-T c cupratesi 1 ' 2 ' 14 

Making use of the solution of Dyson's equation and the results obtained in Appendix A regarding the bare Green's 
function, we numerically determined and plotted A(r, uj) for UDW in Fig. [TJ The four panels show the energy and 
position dependence of the spectral function around the impurity. For example the top left panel corresponds to the 
m = chain, that is the one where the impurity resides. The other three demonstrate the effect of the scatterer on 
the neighboring chains with m = ±1, m — ±2 and m = ±3, respectively. It is immediately clear from the figures 
that the presence of the local perturbation violates particle-hole symmetry. This is a common feature of potential 
scattering and is familiar in dSC 4 ^ and dDW— £ as well. In addition to that, the local electronic structure exhibits 
quite similar patterns on every chain. The fine details and features of these patterns can be nicely separated into 
three distinct components, each having its own microscopic origin. These will be studied in the followings: 

(i) Perhaps the most apparent flavor of the images in Fig. [1] are those curved stripes or waves that are becoming 
ever denser at high energies. They are essentially nothing else but the electronic wavefunctions and can be obtained 
even in a strictly one-dimensional metal. A simple de Broglie picture can already account for the observed periodicity 
along the chain: A = 2n/p, where p — (ui + vpkp)/vF is the momentum of the electron with energy u> measured from 
the chemical potential in a linearized band. The very same electronic waves were found in the spectral function of 
one-dimensional conventional CDWj^ Namely, in Ref. [36[ the effect of open boundary on CDW was studied. Among 
others it was found that the position of zeros in the STM image is determined by kpx = niTVFkF/(u + vpkp), with n 
a natural number. This result is in complete agreement with our simple reasoning based on de Broglie formula. Each 
stripe can therefore be assigned a natural number n (see especially the top left panel with the added curves), though 
they do not necessarily indicate zeros in UDW. This is because the pattern in the present case depends (weakly) on 
the scattering amplitudes too, and an exact agreement is achieved only in the limit U(Q) — U(Q) — > oo corresponding 
to the case of open boundary. 

(ii) An other interesting feature of the STM images is that there is a modulated behavior along the chains with a 
much larger wavelength than A. Namely, A(r, ui) "periodically" takes on the metallic value No, in other words exhibits 
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a beat. This property arises from the position dependence of the zeroth order Green's function, see Appendix A. One 
can check either from the numerics or from the analytical calculation, that on chain m whenever J^ n (2x/£ i ) — 0, the 
LDOS is that of a normal metal. Here J m (z) is the Bessel function of the first kind and £ = vp/h is the characteristic 
lengthscale originating from finite interchain coupling. In quasi-one dimensional density wave materials tb is usually 
of the order of the energy gap A, and this leads to a £ that is comparable to the DW coherence length £o = vp / A, in 
any case it is much larger than atomic distances. This modulated behavior related to the zeros of the Bessel function 
is clearly missing in strictly one-dimensional models and is a speciality of real quasi-one dimensional systems, let it 
be either a normal metal, a CDW or even UDW. 

(iii) The most important property of the electronic structure, at least from the UDW point of view, is definitely 
the low energy behavior around the Fermi level. The subgap structure of LDOS, as shown enlarged in the left panel 
of Fig. [21 is qualitatively different from that of a fully gapped DW^ ' 19 ' 20 i 21 i 36 Namely, in UDW considerable amount 
of spectral weight is accumulated in the intragap regime in the form of two virtually bound quasiparticle states. The 
energies of these impurity states are determined by the poles of the T-matrix 

[i - u ±9l {n ± )f + [u ±92 {n ± )f = o, (12) 

where U± — U(0) ± U(Q). In principle, the solutions of Eq. Q12j) are complex, indicating the resonant nature of the 
virtual states. Making use of the explicit forms of given by Eqs. (|A8|) and (|A9[) . the energies Q! ± and the decay 
rates Vl'± are obtained as 

fi ± EE (4 - iQfi = --A \ p- ( 1 + jl gg^±U , (13) 

± ± ± AW ± ln(4iVo|t/ ± |) V 2 1n(47V |C/±|)y ' V ' 

where we have assumed the impurity scattering to be close enough to the unitary limit so that the result can be 
computed to logarithmic accuracy with ln(4/Vo|C/±|) S> 1. It is only in this limit that the two bound states are 
well defined with f2± <C |J2±|. It is also easy to see, that the finite lifetime results from the finite density of states 
in the subgap coming from nodal quasiparticles (A°(ui) = —(2/7r)g2(u>) ~ \uj\ for small to). In contrast to this, in 
conventional DW with constant gap the binding energies of the corresponding impurity states are purely real leading 
to infinitely sharp resonances and well-defined undamped states. In particular, for example if forward scattering is 
ignored (t/(0) = 0), two symmetrically placed states with opposite energies are foundj 18 ' 21 Equation (fl3)) and the 
presence of impurity induced quasiparticle states in the UDW gap are in precise agreement with that found in d-wave 
superconductor^, and in d-density waveii^ However, in dDW only one such impurity state has been found, which 
is due to the limitation of the strictly point-like impurity potential applied there. Indeed, for such a potential U— = 
and the pole structure of the T-matrix exhibits a single resonance only. Note also, that in dSC and in the PF scenario^ 
of the PG phase, a double peaked resonance is found: one state on both positive (electron) and negative (hole) biases. 
This structure arises from the particle-hole mixing, an essential feature of pairing in SC, and has nothing to do with 
different forward and backscattering. 

In quasi-one dimensional (U)DW systems the fine tuning of chemical potential (doping) does not play such an 
important role as for instance in quasi-two dimensional dDW. It is because it leaves the nesting property unaffected, 
as /i can be incorporated in kp. Consequently, the energies of the intragap resonance peaks are not affected by doping, 
as opposed to that found in dDW, where it scales with /z4 In this respect UDW behaves much the same way as a 
dSG^: there due to the pairing mechanism the quasiparticles are always excited with respect to the Fermi energy, 
leading to impurity induced states with energy pinned almost on the Fermi surface and essentially unaffected by \i. 



IV. FRIEDEL OSCILLATION 



This section is devoted to the analysis of density oscillations caused by the impurity in the Born limit. The 
investigation of this issue is motivated by the findings of Ref. [l8| , valid in one-dimensional conventional charge density 
wave. Namely, it was found, that at zero temperature below the coherence length £o = vf/A, the charge density 
around the impurity is just the sum of the contributions corresponding to the CDW and the Friedel oscillations. 
Beyond £o however, as the necessary electron-hole pairs with energy smaller than 2A are not available, Friedel 
oscillations essentially cannot build up. The exponential decay of the oscillations at this lengthscale was interpreted 
as a tunneling effect. 

Our aim in this section is to perform an analogue but finite temperature calculation in UDW, and find the total 
density below and beyond the coherence length. This enterprise constitutes a more difficult problem than the determi- 
nation of the one-dimensional T = CDW response performed in Ref. [HI because of the following three reasons: (i) 
As UDW is a quasi-one dimensional structure, the finite interchain coupling has to be taken into account as well. This 
introduces an additional lengthscale £ = vp/t^ in the analytical calculation besides £o and the atomic distances vp/D, 
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b. Consequently the Friedel oscillations will become essentially two-dimensional and the chains with index m ^ 
will be affected by the impurity too. (ii) In contrast to fully gapped DW, an UDW is characterized by a momentum 
dependent order parameter and this extra k-dependence makes explicit calculation more cumbersome, (iii) At finite 
temperature T the temperature itself introduces a characteristic length: £1 = vp/T. Although the treatment of finite 
T obviously leads to extra complications in the actual calculation, we try to incorporate its effect in our theory. We 
will see shortly that it is indeed important because it leads to qualitative changes compared to the zero temperature 
results. 

Since in UDW no periodic modulation of either charge or spin density is present, we expect robust Friedel oscillations 
showing up below £q- On the other hand, as UDW is gapless, and nodal excitations are available with arbitrary small 
energy, we expect the oscillations beyond the coherence length to exhibit power law behavior as opposed to exponential 
decay. 

A. Normal metal 

To get started, we first present the results for the density oscillations obtained for a quasi-one dimensional normal 
metal in the Born limit, where Tyg(iw n ) = t 7< j. Using Eq. {8]| and the Matsubara versions of the bare Green's function 
given in Appendix A, for \x\ 3> vp/D one obtains 

„(*, m) = no - ^(QM-l)^ (y ) P (^) ^g^, (14) 

where riQ — kp/(-irbc) is the homogeneous density, P(z) — zsinh _1 (z) and J m (z) is the Bessel function of the first kind. 
This is a very instructive result and it is worth stopping here for a moment and analyze it in detail. The first thing it 
tells us is that in a normal metal in the Born limit only backscattering contributes to Friedel oscillations as the forward 
scattering amplitude drops out from the calculation. It is also clear from Eq. (fT"4|) that the different lengthscales are 
factorized in the sense that they appear in different terms of a product. The last term is the familiar Friedel oscillation 
with periodicity 2kp and algebraic asymptotics (~ M -1 )- This behavior is known to be the consequence of the sharp 
Fermi surface of normal metal at zero temperature. In other words, from a more physical point of view, the long range 
oscillations develop because it is not possible to construct a smooth function out of the restricted set of wave vectors 
|fc| < kp. Here, in our finite temperature result we find that this oscillating term is modulated by a smooth envelope 
P(z), that in the zero temperature limit correctly simplifies to P(z — > 0) = 1. On the other hand, at arbitrary small 
but finite temperature the long-range behavior is replaced by an exponential decay according to P(z 3> 1) ~ 2ze _z .— 
This qualitative change arises from the fact that at finite temperature the Fermi surface is smeared over a thickness 
T in energy and the electronic distribution function becomes smooth and analytic. 

As to the effect of impurity on the parallel chains, it has been taken into account in Eq. (fT4)) by the Bessel function. 
In the extreme limit of decoupled one-dimensional chains, where tf, — > 0, one readily finds J m (2x/t;) — > 8 m Q, indicating 
the fact that screening takes place on the chain only where the impurity resides and all the others are completely 
unaffected. At finite interchain coupling the square of the Bessel function serves as a modulating function: due to its 
quasiperiodic zeros it results in a beat in the induced charge density with wavelength A f=a 7r£, that is certainly much 
larger than that of the Friedel oscillation, n/kp. 

B. Conventional CDW 

The corresponding formula for the impurity induced charge response in a quasi-one dimensional conventional CDW 
can also be obtained in closed form at distances larger than the atomic length scale. Below the transition temperature 
the order parameter is finite, and one can easily show that in this case the Born calculation acquires an extra 
contribution from forward scattering processes. However, it is small and does not exhibit the relevant 2kF periodicity 
characteristic to Friedel oscillation. We therefore omit it here and concentrate on the effect of backscattering only. 
With all this, the analogue of Eq. (|14[) in CDW at zero temperature reads as 

n(x, m) = n - co S (2k F x + <f>) - n A C/(Q)^(-l) m J 2 m (j\ F (^j , (15) 

where 

F(z) = 2K x {z) - 2e* cos(</>) (| - zK (z) - |z (L x (z)K (z) + L {z)K x (z)j) . (16) 
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FIG. 3: The impurity induced electronic density is shown for a conventional CDW (solid) and an UDW (dashed-dotted) on 
the m — chain. It is calculated from Eqs. (|15[1 and (|17[1 , For plotting feF^o = = 10 were applied and n°(x, 0) denotes the 
unperturbed density without impurity. The CDW response freezes out exponentially at large distances. The UDW contribution 
on the other hand is much larger which leads to the fact that the beat property is well observable too, in the present case with 
k F X/2 w 7rfc F £/2 w 15. 



Here n\ = A/|<?| with g being the density wave coupling constant. Furthermore K n (z) and L n (z) are, respectively, 
the modified Bessel function of the second kind and the Struve function. The first two terms together in Eq. (|15p 
is the usual CDW charge modulation. In a pure CDW without impurity, the phase <fi of the order parameter is 
unrestricted due to incommensurability. However, introducing the perturbing potential the phase and therefore the 
overall position of the condensate gets pinned to <j>(U(Q) > 0) = 0, or </>(i7(Q) < 0) = 7r. ls In either case F(z) in 
Eq. (fTfJ)) remains the same, and within the CDW coherence length it simplifies to F(z -C 1) w 2z~ 1 . In this region 
therefore the contribution of the impurity to the Friedel oscillations is precisely that of the normal metal. On the 
other hand at distances much larger than £o we have F(z S> 1) ~ \/2~7rz~ 3 / 2 e~ z , and exponential decay is obtained 
in agreement with Ref . , see also Fig. [3] 

Strictly speaking, Eq. (fT5)) is valid only at absolute zero temperature. However, its validity still holds for low 
temperatures where /3A(T) = £i/£o ^ 1- The justification of this statement is as follows: As we have just seen, 
at large distances exponential decay is obtained with characteristic length £o- This behavior is very similar to the 
effect of finite T in normal metal. The obvious parallelism can be easily understood by remembering the fact that in 
a conventional CDW with constant gap the Fermi surface is smeared over an energy width A. Thus the electronic 
distribution becomes analytic even in the ground state leading to the aforementioned exponential behavior. This 
result is precisely the same as found in BCS superconductor in Ref. [38l |. Raising the temperature introduces a further 
exponential cutoff (~ exp[— 27r|x|/£i]), just like it did for normal metal, that has clearly no observable effect as long 
as £i ^> £o- The opposite relation on the other hand would lead essentially to the normal state. 

Finite interchain coupling tb results in a quasi-one dimensional CDW structure. Its effect coincides exactly with 
that found for normal metal because of the same factors of Bessel functions appearing in Eqs. fTf]) and (|15[) . For 
details see the last paragraph of subsection A. 



C. Unconventional DW 



Now we finally turn our attention to UDW and for the same reasons pointed out in the case of CDW, we again 
consider the effect of the relevant backscattering only. We begin with the zero temperature case. In the ground state a 
closed solution like the ones in Eqs. (fT4|) and (fT5|) valid for all \x\ ^S> vp/D cannot be obtained. Nevertheless, we have 
just seen during the calculation of the CDW response that within the coherence length the metallic result applies. 
This is because at such distances only the high energy electron- hole pair excitations contribute to density oscillations, 
and far from the Fermi energy a fully gapped CDW behaves the same as a normal metal. This latter statement is 
equally true for an UDW as well. Therefore we conclude, that in quasi-one dimensional UDW ground state within the 
coherence length the density oscillations are the same as that of a normal metal given by the zero temperature limit 
of Eq. (fT4|) . Note here, that in UDW due to the vanishing average of the gap over the Fermi surface, the anomalous 
contribution to the total density due to the condensate is missing, n\ = 0. On the other hand, at large distances 
the small energy nodal excitations dominate the static charge response. To make the picture whole, for \x\ 3> £o 
asymptotic expansion to leading order yields 



2 ,'2x 7T ^ cos(2fc F ir) 4z 
n(x,m) — n a = — noJV (7(Q)(— 1) cos — — m 1 



2 



i 2 2fc F £ (2z 2 + m- / 

(17) 

2x 7T \ cos(2A;Fa; + 2(j>) 2m 2 



V ? 2 / 2fc F £ z (2z 2 

where z — \x\/£o- Equation (fT7|) certifies our expectations about the power law decay (~ r~ 3 ) at large distances. At 
this point we would like to call the attention to the following: The argumentation (about the smeared Fermi surface 
leading to exponential decay) we applied in the CDW in subsection B is not applicable directly for an unconventional 
density wave, though a finite order parameter exists here too. In fact the UDW gap is momentum dependent and 
possesses nodes. Around these nodes the electronic distribution function (a^a^) — 2 _1 (1 — £(k) / yj ^(k) 2 + |A(k)| 2 ) 
simplifies to the sharp step function. At other portions of the Fermi surface however, where |A(k)| is close to 
maximum, the distribution function is smeared out much the same way as in normal CDW. Nevertheless, due to 
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FIG. 4: The local density of states is shown at the impurity site versus energy and the impurity strength, where U+ = 
U(0) + U(Q). 

nodal quasiparticles the distribution function remains non-analytical in momentum space leading to the observed 
long-range power law decay. 

At finite temperature charge response is affected only at distances \x\ >> £1, where the factor exp[— 27r|x|/£i] 
becomes dominant. But as long as (3A(T) = £i/£o 3> 1: the power law behavior can in principle be observed in the 
range £o "C \x\ <C £i- In Fig. [3] the density oscillations along the m — chain are compared in a conventional CDW 
and an UDW. Beyond the coherence length the CDW response is hardly observable while it is considerably larger 
in the UDW. Furthermore the aforementioned beat due to finite interchain coupling can be clearly seen in the latter 
case too. This complex behavior signals the presence of the nodal density wave, and seems to be more accessible in 
experiments than the detection of an exponential decay in normal CDW. In the latter case of course, clear CDW 
background is present too, that is certainly measurable by other means. In UDW candidates however, measuring 
Friedel oscillations for example in STM measurements might serve as a useful tool in identifying the low temperature 
phase and to reveal hidden order. 

V. THERMODYNAMICS 

The effect of the single non-magnetic impurity on the thermodynamics can be derived from Eq. ([9]), which gives 
the change in the grand canonical potential in the presence of interaction. This equation will be utilized in this 
section in order to calculate the entropy and specific heat contribution of dilute impurities. According to Eq. |9|). 
we need the exact Green's matrix at the impurity site and at r = —0 imaginary time. The Matsubara summation 
is transformed, as usual by deforming the contour, into a frequency integration along the real axis. This form now 
involves the spectral function at the impurity site, and we get 

G Q/3 (0,0;r<0) = i f dw f{u))e-" T [A+{<j) + apA_ (w)], a,f3 = R(= +1), L{= -1), (18) 

^ J — oo 

where f(ui) is the Fermi function, A + (uj) = A(r = 0, w), and A_ has the same functional form as A+, but instead of 
U+ it involves J7_, that is 

A±{u) = -- - - g + 2g 1 {u ± f\d UJ g 1 (iu ± )\- 1 5{L, u ±). (19) 

U±g 1 y + (U ± g 2 y 



Here t/1.2 stand for the real and imaginary parts of the local zeroth order Green's function given by Eqs. (|A6[) and 
(IA7[) in Appendix A. The LDOS right at the impurity site A(r = 0, u>) is shown in Fig. |4j In addition, ui± denote the 
binding energies of the high energy bound states outside the band determined from 1 — U±gi{oj±) = 0. This is nothing 
else but Eq. (|12| for frequencies where g 2 = 0, i.e. outside the band. It provides us those poles of the T-matrix that are 
located on the real frequency axis. With the aid of Eqs. (|A10[) and (|A11|) the solution reads as cu± = D coth(l / NqII±) . 
These states are indeed very far from the Fermi level, their contribution to the low temperature behavior is obviously 
negligible. Doing so and keeping only the first term in Eq. (fl9|) . we are able to perform the coupling constant integral 
in Eq. ^ analytically, and we end up with 



SQ = -- I dw/(w) 

7T 



arctan [ W±MEtZlL) + arctan ( ^ + ™ U ~ ~ M + 2 arctan ( * 

92 J V 92 J \92 



(20) 



We note at this point that the temperature enters this expression not only through the Fermi function, but also 
through the UDW order parameter A(T) appearing in (?i.2- Nevertheless, the low temperature correction to A(0) 
is small, it is of the order of T 3 ^ and thus cannot affect the expected T linear entropy and specific heati 30 i 31 ' 32 A 
Bethe-Sommerfeld series expansion at low T finally yields the specific heat contribution as 

2^1_AW(0) 

3 V A(0) ' v ' 

where 1/V is the impurity concentration in this single impurity problem. It should be replaced by n.j = Ni/V in the 
extreme dilute limit, where TV, is the number of impurity atoms in the sample. Equation (|2ip shows that the T linear 
contribution is due to forward scattering. Furthermore, niN U (0) / A plays the role of the residual density of states 
at the Fermi level in an impure UDW. These findings are in accordance with Ref. [32| . 
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VI. CONCLUSION 



We have studied the effect of a single non-magnetic impurity in a quasi-one dimensional unconventional density 
wave. In particular, its effect on the local electronic states and on the Friedel oscillations were explored in great 
detail. The potential scattering was treated within the T-matrix approach, with the allowance for different forward 
and backward scattering. In this respect, we extended the widely applied strictly point- like impurity picture. The 
impurity induced local density of states becomes asymmetric with respect to the Fermi energy, signalling the violation 
of particle-hole symmetry. We found a double peaked quasiparticle resonance in the subgap, calculated the energy and 
lifetime of the virtual states, and determined the scanning tunneling microscopy image along the neighboring chains. 
The double peaked nature of the virtual states stems from the generalization for different scattering amplitudes. 
Indeed, for equal forward and backscattering corresponding to a delta potential, we obtain a single impurity state in 
accordance with that found in quasi-two dimensional dDW. The energies and decay rates are found the same as in 
dDW or dSC, which is not surprising at all, as these systems all have a quasiparticle density of states that is linear 
in energy around the Fermi level arising from nodal excitations. The electronic states around the impurity can be 
studied experimentally with STM spectroscopy. This diagnostic tool is able to resolve both the energy and spatial 
dependence of the local density of states by directly measuring the tunneling conductance. To this end, we calculated 
the expected STM image of an UDW as a function of energy and position. 

An external perturbation like an added impurity, or a constraint on the electronic wave functions in the form of 
boundary conditions, are known to cause Friedel oscillations in the charge distribution. In pure UDW no modulation 
of either charge or spin is present. Therefore, robust Friedel oscillations were expected to show up below the density 
wave coherence length. We found indeed, that in this length scale the density oscillations were those of a normal 
metal. This is to be contrasted with the CDW result, where the impurity induced oscillations are superimposed on 
the usual CDW background. On the other hand, contrary to the exponential decay of fully gapped CDW, beyond the 
coherence length power law behavior was expected and found. It is due to the fact, that in UDW nodal excitations 
(electron-hole pair) are available at arbitrarily small energy, thus density oscillations can build up. This algebraic 
behavior at large distances signals the presence of the UDW condensate and could be more accessible in experiments 
than the exponential decay of normal CDW. In addition to that, the possibility of experimental detection is further 
supported by the fact that an UDW does not exhibit static charge density wave background that could overwhelm 
the impurity contribution. 

At last, we calculated the change in the grand canonical potential, the entropy and specific heat contribution of 
the scalar impurity embedded in the UDW host. The calculation was done to infinite order in the interaction. At 
sufficiently low temperature forward scattering produces metallic T-linear behavior, because it breaks particle-hole 
symmetry and causes finite residual density of states at the Fermi energy. 
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APPENDIX A 



In this Appendix we calculate the zeroth order retarded Green's matrix G^(r,r';w) of an unconventional density 
wave in real space. With the aid of the definition in Eq. ([5]), and using the form of left-, and right-moving fields, see 
Eq. ((4]), one obtains 

G%(r, r'; u) = e i(Q/2)[(«-i)r-(0-i)r']+i(^ a -^)^ (r _ ^ ^ a ^ = R{= L{ _ _^ (A1) 

where <fi = (f>R — cj>i, is the phase of the order parameter. As we detached this complex phase factor from the order 
parameter, it is now purely real: A(k) = Asm(bk y ). Furthermore, the homogeneous part reads as 

f<0 ( r . A - 1 V r *kr W + £( k ^ + A ( k )^ (A^ 

^ k ^ 2 -£(k) 2 -A(k) 2 +signHzO' [AZ) 

with pi the Pauli matrices, and the prime on the momentum integration again denotes that the sum is cut off as 
\k x — kp\ < k e , where vpfc c = D is half the bandwidth. We note that in a strictly one-dimensional (no perpendicular 
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coupling between parallel chains, t& — 0) conventional density wave with constant gap, one only needs to perform the 
substitution sin(6fc y ) — > 1 in order to get the desired results. In particular, for the diagonal components one finds 



G° aa (x,u>) c 



1 x ( ik F + 

-e 



A-KVp 



VA 2 



Ex — (D - i^/A 2 -uj 2 ) - Ei ( —(-D - z^A 2 -^) - 2iTrQ(-x)e(D' - U 



1 a: ( zfcp - 

-e 



Va 2 



(A3) 



4ttv f " V \/A 2 - 

Ei ( — (£> + i^A 2 -uj 2 )) - Ei (—(-D + i^A 2 -oj 2 ) ) + 2mQ{x)Q{D' - U) 



and the offdiagonal components read as 

1 iA 



G° a {x,uj) c 



4:1tvf \/A 2 — w 2 

E x ( — (D - 1V /A 2 -^ - Ei ( — (-D - i^jA 2 - uj 2 )) - 2iirQ(-x)e(D / - U) 



1 zA x ( ikw- 

=e 



47rw F VA 2 - 



(A4) 



Ex — (D + iVA 2 -u 2 ) ) _ Ei [ — (-D + zVA 2 -w 2 ) + 2m@(x)Q(D' - U) 

V F I \VF 



Here Ei(z) is the exponential integral, u> means u> + iO under the squareroots, and D' = y/D 2 + A 2 . These formulas 
reproduce those of Ref. [l8| in the limit D — > oo. Now it is convenient to point out, that the corresponding results 
for a strictly one-dimensional normal metal can be obtained easily with the substitution A = 0, and one obtains the 
usual diagonal and translational invariant solution. 

After introducing the strictly one-dimensional results, we now turn our attention to quasi-one dimensional systems 
with finite interchain coupling, % ^ 0. These include the normal metal, a conventional DW, or UDW. We are primarily 
interested in the latter system, because this is the one that exists only in dimensions greater than one. However, 
as the UDW formalism is the most general with its momentum dependent gap, it incorporates the results for each 
case. Quasi-one dimensionality means that there is at least one more (perpendicular) direction in real space, where 
the model is discrete rather than being continuous as it is in the chain direction x. In our notation this additional 
dimension was chosen to be the y direction, where the order parameter varies. In this respect, the position argument 
of the Green's function is r = (x,mb, 0), where b is the corresponding lattice constant and m is an integer indexing 
parallel chains. With all this, the analogous results for the bare Green's function in UDW are obtained from Eqs. (|A3[) 
and (|A4[) by performing first the substitutions A — > Asin(6fcj,), ±D — ► ±D — 24 cos(bk y ), and then integrating over 
the Fermi surface 

d(bk y ) i(6 fc m+ y ca8 ( &fc )) a , , 



(A5) 



where £ = Vp/tb- The integration in Eq. (|A5[) can be performed for a normal metal and a conventional density wave 
in the \x\ 3> vf/D limit, that is for distances much larger than the atomic lengthscale. It leads to the appearance 
of J m (2x/£ t ), the Bessel function of the first kind. For UDW, since the integration cannot be carried out, the Bessel 
function does not show up explicitly. Nevertheless its properties are coded in the integral and we will indeed encounter 
some of them during the calculation of the spectral function and Friedel oscillations. 

We end this Appendix with the presentation of the local Green's function in UDW. At the impurity site x = 0, 
m = the averaging over the Fermi surface can be done analytically and one finds G° a a = S a p(gi + 2.92), where 



ffl M = ReG%(0,0; W ) 



A /D 2 A 2 -w 2 (D 2 +A 2 - 
ZVA 2 -^ 2 



NoF ( £, 



< M < A, 



A < U < D, 



AosignHK (^S) 
^osignHF (-rfLp , 



(A6) 



D < \uj\ < \JD 2 + A 2 , 



VD 2 + A 2 < M, 
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g 2 (Lj)=hnG° RR (0,0;cj) 



-*o*(#) 



-N 




< |w| < A, 
A < \u\ < D, 
D <\oj\< 



VD 2 + A 2 < 



A 2 , 



(A7) 



Here No = l/(TTVFbc) is the normal state density of states at the Fermi level (per unit volume and for one spin 
direction), and K(k) and F(z,k) are, respectively, the complete and incomplete elliptic integrals of the first kind. 
During the calculation of Eqs. (|A6|) and (|A7[) a correction of order 0(tb/D) was neglected as it only affects the 
frequency dependence at high energies around the band edges, and is therefore irrelevant close to the Fermi level. In 
addition to that, in this case the bandwidth can also be taken to infinity and the above formulas simplify to 



.91 



.92 




-iV 9(A - 

-N fe(A - M)M# (£) + e(M - a)k 



(A8) 
(A9) 



On the other hand, if we are interested in the high energy behavior close to band edges, A can be neglected compared 
to D and Eqs. (|A"6j) and |(A"7|) simplify to 



N . 
.9i = — In 

nNo 

.92 = 7T- 



uo + D 



uj-D 



e(i>-H). 



(A10) 
(All) 
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